Spin thermoelectrics in a disordered two-dimensional Fermi gas 
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We study the connection between the spin-heat and spin-charge response in a two-dimensional 
disordered Fermi gas with spin-orbit coupling. We show that the ratio between the above responses 
can be expressed as the thermopower S = — (7rfc_B) 2 Ta'/3ecr times a number R B which depends on 
the strength and type of the spin-orbit couplings considered. We illustrate the general results by 
examining different two-dimensional electron or hole systems with different and competing spin-orbit 
mechanisms and conclude that a metallic system could prove much more efficient as a heat-to-spin 
than as a heat-to-charge converter. 



The moving carriers in a metallic system, electrons or 
holes, transport both electric charge and heat. This gives 
rise to a number of thermoelectric effects as well as a 
deep connection between thermal and electrical conduc- 
tivities. A well known example is the Wiedemann-Franz 
law, which states that the ratio of the thermal to the 
electrical conductivity is the temperature times a univer- 
sal number, the Lorenz number L — ir 2 kg/3e 2 , where ks 
and e are the Boltzmann constant and the unit charge. 
From the theoretical standpoint the validity of the above 
law relies on the single-particle description of transport, 
on the Fermi statistics of carriers and on the assump- 
tion of purely elastic scattering p] [2]. When electron- 
electron interaction is present as in a Fermi liquid, this 
law still holds provided the quasiparticles do not ex- 
change energy during collisions. At low temperatures the 
combination of electron-electron interaction and disorder 
may change this picture [SHU]- Additionally, a magnetic 
field affects both thermal and electrical transport yielding 
both galvanomagnetic and thermomagnetic effects [10] . 
The above situation gets even more complicated when a 
third quantity transported by the carriers - the spin - 
is connected to the previous two by spin-orbit (SO) cou- 
pling. On the bright side, such a connection also opens 
up a plethora of new possibilities related to the manipu- 
lation of the additional spin degrees of freedom. This is 
testified by the recent rapid development of spintronics 
[TT1 [T2"] and spin caloritronics [T3] ■ A fundamental goal 
of spintronics is the achievement of all-electrical control 
of the carriers' spin, made possible by SO coupling as 
exemplified by the spin Hall effect 14-18] . Similarly, an 
important goal of spin caloritronics is the manipulation 
of the spin degrees of freedom via thermal gradients [Hil - 
123] . particularly relevant when energy efficiency issues 
are considered 13J . In this context a noteworthy phe- 
nomenon is the spin Seebeck effect: a spin current ther- 
mally generated in a (metallic or insulating) ferromagnet 
is injected into a normal metal and there, via the inverse 
spin Hall effect, it generates an observable voltage drop 
in the direction orthogonal to the applied thermal gra- 
dient [24H26] . In this much studied case phonons and 



magnons play the leading roles [2TH2TJ] . There are on the 
other hand only few works on thermo-spin transport due 
to the charge carriers' dynamics j2Ql [21], and we wish 
to address this issue considering disordered Fermi gases 
with SO coupling. We will see that a general relation 
between the spin-heat and spin-charge response of such 
systems can be obtained, with the same range of appli- 
cability of the Wiedemann-Franz law. Moreover, we will 
discuss the particular case of the thermo-spin Hall effect 
- the generation of a spin current transverse to a ther- 
mal gradient, also called the spin Nernst effect. In so 
doing we will show that a simple relation connects the 
spin thermopower - the ratio between the spin response 
to a thermal gradient and that to an electric field - to the 
standard electric thermopower, and that the former can 
be strongly enhanced by the interplay between different 
SO coupling mechanisms. 

Let us start with some basic phcnomenological consid- 
erations along the lines of Refs. [3T] , and consider the 
bare-bones situation of an inversion symmetric, homoge- 
neous and non-ferromagnetic material in the absence of 
magnetic fields. A particle current j x can be driven ei- 
ther by an electric field or by a temperature gradient, and 
within the standard semiclassical approach one writes 

m 

j x = L n E x + L 12 (-V X T) = <jE x - e£Ta'(-V x T). (1) 

Here a = —2eNoD is the Drude conductivity up to a 
charge — e, with Nq the density of states at the Fermi 
energy and D the diffusion constant, and a' = d^a, fi 
being the chemical potential. The ratio S = L 12 /Lii is 
the electric thermopower. In the present simple case the 
connection between spin and particle currents due to SO 
coupling reads [3TJ 

f y = -13x = L s n E x + L S 12 (-V X T). (2) 

Here j* is the z-polarized spin current flowing in the y 
direction arising in response to the particle current j x , 
and 7 <C 1 is a dimensionless SO coupling constant. As 
an immediate consequence of Eqs. ([TJ and the spin 



2 



thermopower S s = if 2 /-^!i is equal to 5, since the SO 
coupling constant 7 does not depend on the sources of 
a given particle current. Eq. ([2| breaks down in the ab- 
sence of inversion symmetry, and in order to see how the 
above simple result is modified in a general situation, and 
to study its dependence on competing SO coupling mech- 
anisms, we now move on to a microscopic treatment. The 
latter is actually independent of dimensions (2D or 3D), 
but in order to fix things we consider a disordered 2D 
Fermi gas in the x-y plane described by the Hamiltonian 



H = f- + V(x) + H s , 
2m 



(3) 



with p the 2D momentum and V^x) the impurity poten- 
tial. For the latter we assume the standard white noise 
disorder model with (V(x)V(x')) = (2ttN t)- 1 5(^-x'), 
whereA^o = m/ (2irh 2 ) and r is the elastic scattering time. 
The SO term H so will have different forms in the various 
cases considered below. In the (linear) Rashba case it 
reads 



p x e 2 



(4) 



with a a coupling constant. We assume the metallic 
regime and weak SO coupling conditions, £f ^> h/r, A so . 
Here cf is the Fermi energy in the absence of disorder 
and SO interaction and A so is the SO splitting due to 
H so . From now on h, ks — 1. The a-polarized spin cur- 
rent flowing in the ^-direction due to a generic thermal 
gradient is 



Jk 



X)M,(-Skr), 



(5) 



where N s h is the spin-heat response tensor. Following 
Ref. [7] the latter is given in terms of the imaginary spin 
current-heat current kernel 



[N 



T : 



lim 



[<2sh(i^)]fei 



»n R , Q. R ={1+iO+ 



(6) 

The spin current operator is given by the standard defini- 
tion j% = (l/2){v k , s a }, Vk and s° being the velocity and 
spin operators. Notice that the particle (charge) current 
operator is {—e)jk = (—e)vk- The heat current in the 
Matsubara representation reads 



j'k(p,e n ,e n + fl u ) 



-v/2]ki 



(7) 



with e„ = irT(2n + l),fi„ = 2-rtTv, and € n+u / 2 = e n + 
Q, v /2. The specific form of v k depends on the choice of 
the SO Hamiltonian. For instance in the Rashba case, 
Eq.Q, we have v x _ y — p x ,y/m T aa v ' x . By using the 
Kubo formula the response kernel is given by 

[Q s h\ti (iO,,) ie n+u/2^r \jkGnjlQn+u] , (8) 



where the Matsubara Green functions Q n = Q(p,e n ), 
Qn+v = G(p,( n + ^i/) are matrices in spin space Q n — 
Gn + 2~2 a Gn a<1 - Their analytic properties are determined 
by the sign of the imaginary frequency, therefore when 
performing the momentum integral one obtains a non- 
zero result only if the frequencies e n + VL U and e n have 
opposite signs, which means that e„ is restricted to the 
range —Cl^ < e n < 0. Exploiting that the external fre- 
quency is going to zero one thus has 



p 



Tr [fkQnJlQn 



2tt 



CT sc]fe/ O + ie™), (9) 



where a sc is the static spin-charge (particle) conductivity. 

Although our treatment is general, to illustrate the 
procedure we now consider in detail the Rashba case as 
an example, and we compute the thermo-spin Hall effect, 
i.e. we look for the z-polarized spin current flowing along 
y generated by a thermal gradient along x. The response 
kernel [Q s h] z yx = Q sH reads 



Q sH iiU v ) = T^2J2ie n+u/2 Tr [j*g n j x g n +u] 

£n P 

with Ql = (g n<+ + Q n ,-)/2 and 0° = 
e z ) a (G n ,+ - Qn,-) A whereas 



Gn, 



p- 



T up- 



-sgn(e„) 



(10) 
(p x 

(11) 



2m 2t 
Here p is the chemical potential and the disorder self- 



energy is evaluated in the Born approximation. Eq.(lO) 



only takes into account the so-called bare bubble. Vertex 
corrections [33] IS] will be considered later. According 
to Eq. ([9]) we now have 



2n 



p 



Tr [f y g n j x g n +»] = v sH {}i + ie n ), (12) 



with a sH (p) the static spin-Hall conductivity from 
Ref. [35]. The thermo-spin Hall conductivity therefore 
reads 



N sff T = - lim 



2ttT 



^2 i<W/2 a sH (fi + ie n ) 



(13) 

and expanding in e n yields (the zero order term vanishes) 



N sH k 

1 * I bare 



Tr 2 Tma 2 T 2 

~3 7T ' 



(14) 



To connect this result with that of Ref. 20J, in which 
~N sH T is computed in the clean limit, r — > 00, we rewrite 



Eq. (11) 



Gn,± = 



P 2 ( \ 1 

— T ap + isgn(e n — 
2m \ 2r 



(15) 
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and note that, as long as the temperature is finite, in the 
t — > oo limit the poles are 2ttT away form the real axis. 
Thus the effective replacement l/2r — > irT in Eq.(14| 
yields the clean limit result 



N 



ma 
12ttT' 



(16) 



in agreement with Ref. [2U]. Let us now discuss the ver- 
tex corrections. Taking them into account corresponds to 
sending j* — > J* , j x — > J x and j x — > j£. At the level of 
the Born approximation cither vertex could be renormal- 
ized: the bubble with J* and j£ or that with j* and J x 
are equivalent. Moreover, since we neglect inelastic pro- 
cesses, J x = ie n+ „/ 2 Jx- For the Rashba case it is known 

0, and thus we immediately 



that J x — 0, i.e. 
obtain 



dressed 



o. 



(17) 



However, notice that Eq. (13) holds for any form of the 
SO interaction term H so , no matter whether of intrinsic 
or extrinsic nature. Therefore, once the spin-Hall conduc- 
tivity a sH of a given system is known, its thermo-spin 
Hall conductivity N sH will follow at once. Even more 
generally, from the Matsubara formulation, Eqs. (|6])-(|9]), 
we conclude that the spin-heat response of a disordered, 
SO coupled Fermi gas in the metallic regime is completely 
determined by its spin-charge response. This result holds 
in 2D and 3D, in the presence of arbitrary elastic scatter- 
ing processes, possibly spin-dependent, and beyond the 
Born approximation, i.e. it has the same range of applica- 
bility of the Wiedemann- Franz law discussed in Ref. [7] . 
This is the first main result of our work, which, after a 
Sommcrfcld expansion, can be written in the very simple 
form 



N sh = -eCTa' sc (n). 



(18) 



In other words the spin thermopower S s = —e£Ta' sc /<j sc 
has the same structure as S = —e£Ta'/cr. Whether a 
direct relation between S s and S exists is however not 
obvious and will be our next concern. 

With this in mind, let us now take -ff so to be linear in 
momentum, in which case the spin continuity equations 
assume a particularly simple form. This allows one to 
easily draw a set of more specific conclusions concern- 
ing the thermo-spin response of the 2D Fermi gas, in 
particular regarding the interplay between different SO 
and scattering mechanisms. To be explicit we take once 
more the disordered Rashba model as an example, and 
consider the presence of extrinsic SO mechanisms and 
(white noise) magnetic impurities. That is, we add to 
the Hamiltonian ([3| the terms 



Hcxtr = ~<T X W(X). P) Kn(x) 



B-er£(x-R0, 
(19) 



where Ao is an effective Compton wavelength, 
(y m (x)K 1 (x')) = [3(2 7 r^ r sf )]- 1 (5(x - x') and r sf 
is the spin-flip time [36] [37] . The s y continuity equation 
reads 



d t s y + V ■ j» = -2ma€ - 



3t s - 



1 

TEY 



(20) 



with tey = t(\qpf /2)~ 4 the Elliot-Yafet spin-relaxation 
time due to H extr - Assuming a homogeneous electric field 
applied in the x-direction, the spin current j* in the dif- 
fusive regime is given by 



jy = 2maDs y - ^aE x , 



(21) 



where 7 = 7i ntr + 7 S j + 7 SS is the SO coupling con- 
stant due to intrinsic and extrinsic mechanisms, with 
7intr = —ma 2 T, 7 S j = (Ao/2) 2 to/t the side-jump contri- 
bution and 7 SS = (Aop_F/4) 2 (27rfV i'o) the skew-scattering 
one. vo is the scattering amplitude (see Ref. [38l [39] 
for details). In a homogeneous bulk in steady state 
the spin-Hall conductivity is easily computed, a sH = 
[1/(1+C)]7<7, where C = t s /t dp , with 1/t dp = (2ma) 2 D 
the Dyakonov-Perel spin-relaxation rate, and 1/t s = 
4/(3r s f) + 1/tey- Via Eq. (18) one concludes 



JsH 



S s 



7_ 
7 



c 

i + C 

7 



sH 



i + C 



(22) 



(23) 



with the spin Hall thermopower S s = N sH /a sH . In the 
above, primed quantities are derivatives with respect to 
the chemical potential /i. Notice that the simple phe- 
nomenological argument of the introduction overlooks 
the /i-dependency of 7: the conclusion S s — S holds 
only for an energy-independent 7. Both a sIi and N sH 
depend on the ratio between td p and r s and are in prin- 
ciple tunable, either by varying the doping, which af- 
fects r s , or by modulating a by varying the gate po- 
tential. Let us consider some interesting cases using 



Eq. (22) and Eq. (23). When only Rashba SO and mag- 



netic impurities arc present, we have t s = 3r s f/4 and 
7 = 7int- By evaluating the various derivatives we ob- 
tain 7' = 0, Q' = C//i, cr' — cr/fi, which gives us the spin 
thermopower 



5, = ~eCT- 



1 



,i+c- (24) 

When SO from impurities is present, too, the terms 



y/7, C'/C m Eq. (22) are modified, leading to 



5, = —e£T- 



1 



Tss 

7 



c 



i + C 



1 - 



2r s 
tey 



(25) 



The results so far obtained can be generalized to include 
the effects of the linear-in-momentum Dresselhaus SO 
term described by the Hamiltonian 



(p x a x - p y a v ) . 



(26) 
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FIG. 1. The spin thermopower S s of a disordered 2D-electron 
gas with numerous competing SO mechanisms. Typical val- 
ues for GaAs quantum wells are: mobility fj, — 10 4 cm 2 /Vs, 
density n = 10 12 cm -2 , effective extrinsic wavelength Ao = 
4.7x 10 _8 cm, Dresselhaus coupling constant ft/3 = 10 _12 eVm. 
There follows j ss » 7mtr,7aj, r E y » t§ p . The Rashba 
coupling constant can be modulated by the gate potential 
[401141) . Each panel shows the ratio S s / S as a function of the 
ratio q//3 for a given Elliot- Yafet scattering strength, strong 
to weak from top left to bottom right - panel 3 corresponds 
to standard GaAs. Magnetic scattering is strongest for the 
dotted curve, r s f /tEp = 1, and strong (weak) for the dashed 
(solid) curves, T sf /rg P = 2,3(10,20,30). 



-mr(a 2 — 

/3 2 ), 1/tdp = (2m) 2 (a 2 + f3 2 ) D = l/rg P + l/rg P and 



It suffices to replace in the above 7i ntr 



r s _ ^ T s I ( T DP T Dp) 

t~dp t s /tup + 1 



(27) 



Derivatives are trivial, but yield expressions too cumber- 
some to be conveniently written down. The results are 
thus plotted in Fig. [TJ and show the sensitivity of the 
spin thermopower to the various physical parameters in 
play. A modest modulation of the Rashba coupling con- 
stant could substantially modify S s , either enhancing or 
decreasing it depending on the systems characteristics - 
we considered ratios a/j3 well within current experimen- 
tal capabilities [40l [41] . We will come back to this point 
in a moment. Let us now consider our final example, a 
2D hole gas as analyzed in Ref. [32]. The SO interaction 
is cubic in momentum 

H so = a H a x [p y (3p 2 - p 2 )] + a H a y [p x (3p 2 - p 2 )] , 

(28) 



and the spin Hall conductivity reads [32] 

V(V i) ± 

(4ry 2 + l) 2 M T 
with 77 = otHPpT [43]. Proceeding as before one gets 



S, = —e£T- 



3 (I2r) 2 - 1) 
(4r? 2 + 1) (Arj 2 - 1) 



All previous result can be cast in the simple form 



(30) 



(31) 



with R s a number which depends on the various compet- 
ing SO mechanisms. Eq. (31 1, which is our second main 
result, looks physically quite reasonable: in a metallic 
system in which electrons (or holes) are the sole carriers 
of charge, spin and heat, the heat-to-spin and heat-to- 
charge (particle) conversions are expected to be closely 
related. The examples considered show however that 
R s > 1 can be easily achieved: in standard GaAs sam- 
ples with Rashba SO and extrinsic mechanisms one may 
estimate R so ~ 3 [55], and the same value is obtained 
in a two-dimensional hole gas with purely cubic Rashba 
SO in the diffusive regime (77 <C 1). If Dresselhaus SO is 
also taken into account, similar values could be achieved, 
as shown in Fig. [T] This suggests that metallic systems, 
typically characterized by low thermoelectric efficiencies, 
could be much more efficient in heat-to-spin conversion 
and therefore play a front role in spin caloritronics. Of 
course, whether substantially higher R s values can be 
reached in different systems, e.g. in transition metals 
which already show a giant spin Hall response [UJ |3S] , 
or more exotic ones such as p-doped graphene [35] or 
topological insulators like HgTe [37] , is an open and rele- 
vant question. Indeed, it would be interesting to establish 
whether it is always possible, within the regime in which 



the general expression (18) holds, to find such a simple 
connection between S s and S. We therefore believe it 
desirable to experimentally test Eq. (31 1. 



In conclusion, we have studied coupled spin and ther- 
mal transport in a disordered and SO coupled Fermi gas, 
and shown the existence of a general expression for the 
spin thermopower S s with the same structure and an 
identical range of validity of Mott's formula for the elec- 
tric thermopower S. Finally, we have derived a simple 
and physically transparent relation connecting the two 
quantities which could be experimentally tested and sug- 
gests that metallic systems could be much more efficient 
in heat-to-spin than in heat-to-charge conversion. 
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